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Abstract 

A revised Marko-Siggia elastic model for DNA double helix [Macromolecules 
27, 981 (1994)] is proposed, which includes the WLC bending energy and 
a new chiral twisting energy term. It is predicted that the mean helical re- 
peat length (HRL) for short DNA rings increases with the decreasing of chain 
length; while for very long chains, their mean HRL is the same, indepen- 
dent of both the chain length and whether the ends are closed, it is longer 
than the value for rectilinear DNAs. Our results are in good agreement with 
experiments. 
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Single- molecule extension experiments on DNA molecules show that freely fluctuating 
open chains (FFOCs) could be well described by the inextensible worm-like chain (WLC) 
model [1-3]. For a WLC chain of total length L, its intrinsic elastic energy is of the form 
(3E = / L AK 2 ds/2, where A ~ 150 base pair (bp) is called the bending persistence length 
and k = \d s t\ is the curvature, the change rate of the tangent unit vector t(s) at arc- length s, 
(3 = 1/ksT with ks being the Boltzmann constant and T the environment temperature [3]. 
However, in addition to bending ones, double stranded DNA molecules have also twisting 
degrees of freedom, and the total intrinsic energy formula for a deformed DNA chain is still 
unclear [4-10]. To know the exact form of the energy formula is very important for the study 
of DNA configurational properties, especially in the case of torsionally constrained DNAs, 
such as covalently closed DNA rings in cells. For example, it has been widely accepted that 
bending elasticity and twisting elasticity determine to a large extent the particular tertiary 
structures of DNA rings [4-6] . Previous studies often regard a DNA chain as a thin elastic 
rod with isotropic cross section [11], the total elastic energy is assumed to be 



with bending and twisting deformations being independent of each other, here C is called the 
twisting persistence length and u is the spatial angular frequency of the unstressed DNA 
double helix [4-6], and Q 2 + Q 2 , = k 2 [11]. Although this simple achiral model is useful in 
some cases and it seems to be the most natural extension of the already verified WLC model, 
it can not properly describe the chiral characteristics of real DNA chains, this chirality of 
DNA molecules was clearly demonstrated by the single-molecule experiment of Strick et al. 
[7]. Another very important phenomenon related to the chirality of DNA molecules is that 
Nature prefers "linking number deficit" in circular DNAs [12-14], this bias can not be well 
explained by the achiral model (1) studied previously [4-6]. 

Recently, a chiral elastic theory was proposed by Marko and Siggia to incorporate cou- 
pling between bending and twisting deformations in the energy formula [8]. This model was 
based on a careful consideration of the intrinsic symmetry of DNA chains. Later on, Kamien 
et al. extended the Marko-Siggia (MS) model to investigate twist— stretch coupling of highly 
extended DNA supercoils and found good agreement with experiment [9, 10]. However, when 
applying the MS model to the case of FFOCs, we find that it is in general not compatible 
with the already verified WLC theory. To overcome this shortcoming, in our present work 
we first propose a revised version of the MS model to ensure this compatibility. It is shown 
that this can be attained if we further hypothesize that unstressed circular DNAs form flat 
circles. The corresponding internal deformation energy is also consisted of two parts as in 
Eq. (1), the bending energy and the twisting energy. The only difference is that Q3 in (1) 
is replaced by f2 3 + (B/C)VL\ in the twisting energy. After deriving the new elastic energy, 
we use this model to discuss the mean helical repeat length (HRL) of open and closed DNA 
chains. Our result show that for short ring-shaped DNAs, the shorter the chain, the longer 
its HRL. This tendency is consistent with experimental observations. For very long chains 
with twisting freedom, i. e., open chains or closed chains with at least one defect, their mean 
HRL is independent of both the chain length and whether the ends are closed or not, its 
value is longer than that of rectilinear DNAs. 

Experimentally, it was discovered that random solution DNAs have a mean HRL signifi- 
cantly longer that of rectilinear DNA fibers [15, 16]. However, the reason for this discrepancy 
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has been obscure for quite a long time. Some researchers suggested that maybe ionic con- 
centrations differ in fibriform and solution DNAs, causing a observable effect on the HRL. 
Our present theoretical results indicate another possibility: that the chirality of DNA might 
be the real reason for this discrepancy. In fact, we find that the ionic conditions in DNA 
fibers with high humidity and in solution DNAs differ only slightly in the experiment of 
Refs. [16] and [15], so we feel our present explanation may be more reasonable. 

First we briefly review the main points of the MS model [8-10]. The configurations of an 
inextensible polymer are specified by three orthonormal unit vectors {u(s), n(s), t(s)} along 
the chain, where t is the axial direction vector of the DNA double helix and u is a unit 
vector perpendicular to t and pointing from one back-bone chain to the other, n = t x u. 
It proves to be convenient to use Euler angles by setting u = ei, n = e2 and t = e3, with 
d s e,i = O x e, (i = 1,2,3), here Q, = (^1,^2,^3) is angular velocity of the frame {e^} [11]. 
Symmetry analysis shows that the polymer's properties should remain unchanged under the 
transformation {ei — > — ei, e 3 — > — e 3 }, thus, under the fundamental assumption that an 
undistorted open DNA is a linear double helix with spatial frequency u , the most general 
elastic energy up to quadratic order in the deformations should be of the following form [8] 

(3E C = j ds[^nl + ^n 2 2 + j(n 3 - ^ ) 2 + - u )] (2) 

the first two terms are related to bending deformations and A', A are bending persistence 
lengths along the directions e 2 and ei, respectively, the third term is twisting energy, the 
last term is caused by bend— twist coupling and B is the coupling constant [8]. Eq. (2) can 
be further extended to include stretch— twist coupling [9, 10], but this effect is not important 
for DNAs at average conditions. 

It is easy to know that, for FFOCs who can twist freely so we need only to consider 
bending deformations, the MS model (2) is not equivalent with the WLC internal energy 
/ A(Vl\ + il 2 ,)ds/2, unless the condition B 2 = (A' - A)C is satisfied perfectly [see Eq. (10)]. 
In other words, if we assume the correctness of model (2) and regard it as the starting 
point, it will not lead to the WLC energy for freely-twisting chains, unless B 2 = (A' — A)C. 
However, the WLC model for open DNA chains has been confirmed by experiments [1-3], 
so for the MS model to be appropriate the bend— twist coupling constant should satisfy this 
strict requirement. However, one would wonder why should the twist— bend coupling B be 
completely determined by bending and twisting characteristic lengths? Is there any intrinsic 
reason? In the following we will see this relation is actually implied in the MS model, but we 
need to add in the MS theory another fundamental assumption that an undistorted closed 
DNA ring should form a fiat circle. 

To show this, let's consider the case of short chains. For a short DNA chain (less than 
about one bending persistence length), entropic elasticity caused by thermal fluctuations is 
neglectable and the chain's free energy equals its internal deformation energy Eq. (2). The 
first variation 5^{(3E C + A ■ J L tds} = with respect to the three Euler angles 9, <fi and 
ip [11] gives the equilibrium shape equation, here A = (Ai,A2,A3) is a Lagrange multiplier 
to take account the possible constraint of fixed end-to-end distance. The shape equation is 
composed of the following three equations [17] 

A' {(f) smipsmO + 9' cos^)0' sin ip cos 9 — A'[((f)' sinipsm9 + 9' cosip) cosip}' 
+A((j)' cos ip sin 9 — 9' sin ip)cf)' cos ip cos 9 + A[(<f>' cos ip sin 9 — 9' simp) sin ■?/>]' 
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-C(<f>' cos 9 + ip' - uj )(f)' sin 6 + B(<f>' cos 9 + ip' - uj )<P' sin ip cos 9 

—B[(4>' cos 9 + ip' — lu ) cos^]' — B(4>' sin ip sin 9 + cosip)(p' sin 9 + Ai sin cos 

— A 2 cos0cos# — A 3 sin# = 0, (3) 



—A' [((/>' sin ip sin + 9' cos ip) sin ^ sin 0]' — A [((/>' cos^ sin 6* — 9' sin^) cos ip sin 9]' 

-C[(4>' cos 9 + ip' - co ) cos 9]' - B[(<f>' cos 9 + ip'- uj ) sin i\) sin 0]' 

—B[((f>' smipsin9 + 6>' cos^) cos#]' + Ai sin 9 cos <p + A2 sin 9 sin <p = 0, (4) 

A'(<f>' sin ^ sin 9 + ^' cos-?/;)(0 / cos-0 sin# — ^'sin^) 
+A(<f>' cos ip sin 9 — 9' sin ip)(—(p' sin ■?/> sin 9 — 9' cos ip) 
-C((p' cos 9 + ip' -u )' - B((p' sin ip sin 9 + 0' cos 

+ J B(0 / cos^ + ^ / -cj o )(0 / cos^sin^-^sin^) = 0. (5) 

In Eqs. (3-5), differentiation with respect to arc length s are denoted by superscript prime (it 
should not be confused with that of the elastic constant A'). Eqs. (3-5) determine bending 
and twisting manners of short DNA chains in equilibrium. To check the validity of these 
equations, we see that (3-5) require that the linear equilibrium shape with 9' = <p' = 
should have ip' = uq, and its HRL is ho = 2ix/ujq. This is in agreement with the previously 
mentioned fundamental assumption of the MS model [8], so Eqs. (3-5) are valid. 

Now we investigate whether a flat-circular configuration can be a equilibrium one. For a 
flat circle (its two Euler angles can be set to be 9 = tt/2 and <p' = p = 2tt/L for convenience) 
to be a equilibrium shape, Eqs. (3-5) require that the third angle, ip, should satisfy the 
following equations, 

(A - A')p(sin ip cos ip)' - Cp(ip' - uj ) - B[{ip' - u ) cosip]' - Bp 2 simp - A 3 = 0, (6) 
(A' — A)p 2 sin ip cos ip — Cip" — B&op cos ip = 0, (7) 
2(A — A ') pip' sin ip cosip — B{ip') 2 cosip + Buj^ip' cosip — Bip" smip 

+Ai cosps + A 2 ships = 0. (8) 

These equations are mutually compatible if and only if [20] 

B 2 = (A' - A)C (9) 

and the Lagrange multiplier A is set to zero. Under condition (9) we can get the sole solution 
of ip to be [17] 

. 2nB , . . 

~ ~CL Sm ^ 

with p in Eqs. (6-8) being replaced by 2ir / L. Thus we know that the necessary and sufficient 
condition for an undistorted closed DNA to form a flat circle is Eq. (9). Note that (9) is 
just also the compatibility condition between the MS model and the WLC model. So 
for the MS model to be appropriate, in addition to the assumption that unstressed linear 
DNA configurations have spatial twisting frequency iv , another fundamental assumption is 
needed. This assumption is that unstressed closed DNA rings form flat circles, under which 
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(9) will certainly be hold. Stable DNA circles have been observed in various experiments 
[18-20], so this assumption is reasonable. 

Correspondingly, we derive the final form of the elastic energy after inserting Eq. (9) into 
Eq. (2), 

(3E C = J ds[^(nl + n 2 2 ) + ^(n 3 -cu + ^) 2 ] (11) 

Eq. (11) is the central result of this paper. It is interesting to see that, compared with the 
previously mentioned achiral model (1), the only new thing of Eq. (11) is that a new term 
(B/C)Qi is added into the twisting energy. In what follows we will discuss the possible 
influences of this new term to DNA twisting manners. 

For short equilibrium DNA rings, from Eq. (10) we know that their HRL is 

/^ = {(^) 2 -(^) 2 }- 1/2 , (12) 

which are longer than ho = 2n/uo, the value for undistorted linear chains. Especially 
interesting of Eq. (12) is that it predicts that the shorter the chain, the longer its HRL. 
Such a tendency in HRL has been observed in various experiments [18, 19], and its real 
nature has been controversial. The present work suggests it to be induced by the twisting 
energy in (11). 

The case of long chains is more important and interesting, but much more difficult to 
tackle. Here we should use statistical methods. Model (11) shows that for a chain with 
twisting freedom its axial (t) distribution is just that of a worm-like chain, with p(t, t , s) oc 
J t * exp[— (A/2) J Q L t' 2 ds], and for each specific axial configuration, (ip')twist = 0'(s) cos 9— 
(B/C)k(s) sin(-0 + a(s)), where (• • •) twist means average respect to twisting distribution and 
a(s) = arctan(#'/0' sin 9). Taking into account the fact that ip changes much faster than 
and 9, we can take neglect the term 0' cos 9 and take k(s) = (0 /2 sin 2 9 + 9' 2 ) 1 ^ 2 as constant 
while calculating ip(s). Thus, (ip) twist — — {B/C)k s sm(ip + const) and the instantaneous 
HRL at arc- length s is h(s) = (2n/u )(l - B 2 k(s) 2 /C 2 u%)~ 1/2 . Consequently the mean HRL 
is 

B 2 

h = (H S ))bend = h [l + 2C 2 U 2 ( K (. S f) bend] (13) 

where (• • -)bend means average with respect to the WLC chain. We see from Eq. (13) that 
h > ho whenever B ^ 0, in qualitative agreement with experiment of Wang [15]. However, 
(n(s) 2 )bend is difficult to calculate for a WLC chain, due to the fact that |t| 2 = 1. Here, we 
have to adapt a self-consistent field method to convert this local constraint to a global one 
such that J L t 2 ds = L and determine the corresponding Lagrange multiplier self-consistently 
[21]. The self-consistent field internal energy is f3E' = J L (A/2)t' 2 + 'jt 2 + A • t]ds where A is 
for possible end constraints. Detailed calculation shows that 

«s) 2 ) bend = ^ (14) 

which is independent of the value of A, i. e., independent of whether the chain is open 
(A = 0) or closed (A ^ 0); and it is independent of chain length as it should be. Thus 
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Eq. (13) predicts that solution DNA, whether linear or circular, has the same mean HRL 
longer than that of rectilinear DNA. 

For very long closed DNA chains (about several ten thousand base pairs), gel elec- 
trophoresis experiments did reveal a significant increase in mean HRL, with h = 10.4 bp 
[15], however, the HRL for rectilinear DNAs is only 10 bp [16]. As mentioned before, for a 
long time the real reason for this phenomenon is not clear. Our present theory can give a 
natural and reasonable explanation, proposing that this discrepancy is induced by the chiral 
twisting energy in model (11). To be more quantitative, we insert the experimental values 
into Eq. (13) and estimate that B/C ~ 17.8. This relatively large value indicates that A' 
is much higher than A in Eq. (2). We wish future experiments will check the results of our 
present theoretical work. 

As mentioned before, another interesting phenomenon observed in experiments is that 
in prokaryotic cells and in some yeast cells many circular plasmid DNA molecules of several 
thousand bp are found with linking numbers typically about 5 per cent less than that of the 
relaxed ring [19]. In fact, almost all naturally occurred DNAs are negatively supercoiled, 
i. e., with deficient linking numbers. Why should natural DNA rings prefer a deficit rather 
than an equilibrium or an excess in linking number? The real nature is still unknown, maybe 
this kind of bias is also related with the chiral characteristics of DNAs described in model 
(11). Investigations based on model (11) with topological constraint in now in progress by 
the present authors. 

In summary, we have proposed a revised Marko-Siggia chiral elastic model for DNA 
molecules and discussed its predictions on DNA double helix mean helical repeat length. 
The theoretical results show that for short DNA rings, their mean HRL increases with the 
decreasing of chain length; while for very long chains, whether open or closed, their mean 
HRL is independent of chain length and is longer the value for rectilinear DNAs. These 
results are in good agreement with experiments. We expect future experiments can check 
our theoretical results. 

This work is supported by the National Natural Science Foundation of China. We would 
like to thank Dr. Liu Quanhui, Dr. Yan Jie and Zhao Wei for helpful communications. 
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